Time-reversal symmetry and random polynomials 



Daniel Braun 1 , Marek Kus 1 ' 2 and Karol Zyczkowski 1,3 
Fachbereich Physik, Universitdt-GH Essen, 45117 Essen, Germany 
2 Center for Theoretical Physics, Polish Academy of Sciences, Warsaw, Poland 
3 Instytut Fizyki, Uniwersytet Jagiellonski, Krakow, Poland 



> 
o 



> 

o 
o 



On 

"c 

i 

o 
o 



- 1—^ 

x 



We analyze the density of roots of random polynomials 
where each complex coefficient is constructed of a random 
modulus and a fixed, deterministic phase. The density of 
roots is shown to possess a singular component only in the 
case for which the phases increase linearly with the index 
of coefficients. This means that, contrary to earlier belief, 
eigenvectors of a typical quantum chaotic system with some 
antiunitary symmetry will not display a clustering curve in 
the stellar representation. Moreover, a class of time-reverse 
invariant quantum systems is shown, for which spectra dis- 
play fluctuations characteristic of orthogonal ensemble, while 
eigenvectors confer to predictions of unitary ensemble. 

The distribution of roots of polynomials of high degree 
with random coefficients was investigated recently in con- 
nection with properties of quantum chaotic systems . 
In particular the authors of the cited references consid- 
ered the coherent state representation of eigenstates of a 
quantum mechanical spin system with the total spin S. 
The polynomials in question have the form 



P(z) = 




N = 2S 



(1) 



where stand for binomial coefficients and at are com- 
ponents of an eigenvector. The complex variable z is 
connected to the Bloch sphere angular variables 9, (f> via 
z = tan(0/2) exp(i^). It was shown by Leboeuf and Voros 
|jj , that for large values of S when the quantum system 
in question is chaotic the distribution of the roots is given 
by 



P{z) 



N 
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corresponding to the uniform distribution of the roots 
over the Bloch sphere. This is the consequence of the fact, 
that in the semiclassical limit N — > oo the components 
with respect to a "generic basis" of the eigenvectors of a 
chaotic systems are independently normally distributed 
(see H and references therein). 

The details of the distribution of the components 
depend on symmetries of the system in question. For 
systems which are not time-reversal invariant the eigen- 
vector components are complex, with independently, nor- 
mally distributed real and imaginary parts, whereas for 
time-reversal invariant systems the eigenvectors can be 
made real (also with normally distributed components). 
In the latter situation the uniform distribution (j^) is 
modified. In particular, the roots tend to concentrate 



on the real line Imz = 0, which is a symmetry line for 
the roots (if zq is a root then it's complex conjugate Zq is 
also a root) see also below). When projected back 

on the sphere the symmetry line is the great circle cj) — 0. 

This simplest situation corresponds to the case when 
the time-reversal operator is represented by the com- 
plex conjugation operator. On the other hand it is 
known that generalized time-reversal symmetries, rep- 
resented by the complex conjugation supplemented by 
a unitary transformation, influence statistical proper- 
ties of eigenvector components in the same way as 
the conventional time-reversal symmetry ||. As an 
illustration the authors of Refs. considered var- 

ious models of the so called kicked top system [||, 
which is described by the one-step evolution opera- 
tor of the form U = exp(— if\) exp(— if 2) exp(— i/3) 
with fi = fi(S x ,S y ,S z ),i = 1,2,3 polynomial func- 
tions of the components of the spin operator S = 
(S x , S y , S z ). The simplest case displaying chaotic dy- 
namics in the classical limit is obtained by choosing 
Uq = exp(—ifiS x ) exp(— ipSl/2S) with appropriate val- 
ues of the parameters /x and p. It has two general- 
ized time-reversal symmetries T\ = exp(—i[j,S x )K and 
T2 = exp(— ii-iS x ) exp(iTrS y ) exp(iirS z )K , TiUT~ x = 
both being compositions of linear rotations with the com- 
plex conjugation operator K. The rotations shift the 
symmetry line from the great circle <j) = to other ones, 
the phenomenon exhibited by the numerical investiga- 
tions performed by the authors of Refs. ^||. 

A non-homogeneous distribution of zeros of Husimi 
functions is linked to statistical properties of coherent 
states expanded in the eigenbasis of the Floquet oper- 
ator. In particular, the number of relevant eigenstates 
and the entropy of coherent states || was found for 
this model to be smaller than average along the symme- 
try lines Ti. A smaller number of significantly occupied 
eigenstates denotes a larger number of weakly occupied 
states, in consistency with investigated clustering of zeros 
of eigenstates in Husimi representation along the symme- 
try curves. Moreover, the distribution of expansion coef- 
ficients of a coherent state localized sufficiently far away 
from the symmetry lines is statistically indistinguishable 
from properties of a generic coherent state of a system 
without any antiunitary symmetry [g] . This corresponds 
to recent result of Prosen j|] , who showed that the den- 
sities of zeros of random polynomials with real and com- 
plex coefficient are equal sufficiently far away from the 
real axis. 
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In order to break the generalized time-reversal symme- 
try the original model Uq was supplemented by an nonlin- 
ear rotation h = qS%/2S (in Ref. §) or f{ = qS 2 z /2S (in 
Ref. H) instead of f\ = Q. In their numerical investiga- 
tions Bogomolny et al. observed vanishing of the concen- 
tration of the roots which they attributed to the breaking 
of the time-reversal symmetry. In what follows we will ar- 
gue that the concentration of the roots on the symmetry 
lines happens in the case of generalized time invariance 
only exceptionally and as such cannot be treated as a 
criterion discriminating between the time reversal invari- 
ant and noninvariant systems. In particular the kicked 
tops Ui = exp(—iqSy)Uo and C/ 2 = exp(—iqS 2 )Uo differ 
with respect to the statistical properties of the spectra 
for generic values of the parameter q. The additional 
rotation term brakes all generalized time-reversal sym- 
metries for the first top and U\ pertains to the circular 
unitary ensemble (CUE), while the second still possess 
such a symmetry 

T' = exp(—iqS 2 ) exp(—ipS x ) exp(—iqSl)K, (3) 

and its spectrum is typical to circular orthogonal ensem- 
ble (COE), irrespective of the value of q. Note that the 
above operator is constructed of a nonlinear unitary ro- 
tation (quadratic term S 2 in the exponent), in contrast 
to the operators T\ and T 2 . 

Inasmuch as level statistics reveals directly the sym- 
metry properties of quantum systems, special care has to 
be taken interpreting the statistical properties of eigen- 
vectors, since their distribution depends on the basis cho- 
sen. For example, the distribution of eigenvectors of Uq 
in the S z basis does not confer to COE predictions. The 
agreement with random matrices is recovered in S x ba- 
sis: the geometric symmetry of the top manifests itself in 
the structure of operator Uq. It splits into two parities of 
size S and 5 + 1, which have to be treated separately to 
achieve results according to random matrices. In earlier 
references [^p|-^2| the variables of the top were exchanged 
x «-» z, what gives the same effect. 

The distribution of eigenvectors can be characterized 
by their mean entropy H, which for random matrices of 
size N is equal to H(N, (3) = ^(N(3/2 + 1) - *(/?/2 + 1), 
where VP stands for the digamma function and (3 — 1 
for COE and (3 = 2 for CUE |]|. Figure 1 presents the 
entropy of eigenvectors relative to the entropy of CUE for 
two tops Ui and U2 as a function of the control parameter 
q. Observe similar behaviour for "unitary" top U\ and 
the "orthogonal" top C/ 2 ! The dips in the data for unitary 
top at q — and q = p = 6.0 correspond to transitions 
to the orthogonal class, while U2 pertains to COE for 
any value of q due to the symmetry (^). This difference 
is visualized in level spacing distribution P(s) displayed 
in the inset. An explanation of this fact is simple: out 
of any "orthogonal" spectrum D\ by a generic unitary 
rotation W one can produce an operator Uw — WD\W^ 



which enjoys COE-like properties of the spectrum and 
CUE- like properties of the eigenvectors. This is exactly 
the case of the top t/ 2 , for which the operator exp(— iqS 2 ) 
plays the role of W . Observe that U2 is similar to the 
orthogonal top C/ 2 — exp(—ip,S x ) exp[— i(p + q)S 2 /2S}. 

A similar effect is visible in the distribution of zeros 
of Husimi function representing eigenvectors: both tops 
show lack of roots concentration lines as shown for U\ in 
H and for U2 in || , even though they belong to different 
universality classes. 

In order to understand the above announced results let 
us derive the density of roots p of a polynomial (Q), where 
ak are Gaussian distributed random quantities with fixed 
but arbitrary phases <fk'- 



a-k = r k e lVk 



the Tk being distributed according to 



P(r k ) = 



1 
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(4) 
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We will make use of the same technique employed in 
Ref. U, namely representing p(r, ip) by the Kac-formula, 



p(z) = S[P(z)} 



dP(z) 



dz 



(6) 



and then expressing the delta-functions for the real - 
and imaginary parts of P(z) as Fourier integrals. We 
then get, in full analogy with Eq.(C6) in Ref. S: 
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klyj 'c k N C l N r k ri r ^+l-2 e i{ Vk - Vl ) e i{k-l) v j . 
exp ^2 r n («n cos (p n + f3 n sin cp n } (7) 



where 



k=0 

On = ir 11 \ZC%(cos(ntp)£i + sin(n<^)^ 2 ) 
f3 n = ir n y/U^(cos(nip)^2 - sin(nv?)£i) 



(8) 
(9) 



and z — re lip . Averaging over the random coefficients 
amounts now to simple Gaussian integrations. The re- 
sulting average density can be cast in the following form: 



(p(r, <p)) = J #1 / d &(A + Sei6 + Cg + D^l) 
exp(-a& - b£ - 2c£i6) , (10) 

where 
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B = 



C = 



D 



Y^k 2 r 2{k - 1] C k N 

k=0 
N 

- h a sin(ip k + kip + tpi + lip) 

k,l 
N 

- h a cos(ip k + kip) cos(ipi + lip) 

k,l 
N 

- hki sm(ip k + kip) sin(<pi + lip) , 



k.l 



and h kl = 



(11) 
(12) 

(13) 
(14) 
-lip). 



fclr^+i-iyCftCkcosfat + hp-ipi 
The coefficients of the quadratic form in the exponential 
are given by 



b = 



\^C^r 2n cos 2 {ip n +nip) 

1 N 

-^C^r 2 "sin 2 (^+^) 

n=0 
1 N 

-^Qr 2 " S in(2(^+n^)). 



n=0 

It has the eigenvalues 



Ai, 2 = {a + b± y/(a - b) 2 + 4c 2 )/2 . 



(15) 
(16) 
(17) 

(18) 



and can be diagonalized by a rotation in the (£i, £2) plane 
by an angle 7 = — arctan((a — b — y(a — fo) 2 + 4c 2 )/2c). 
The ^-integrals are Gaussian again and lead to the fol- 
lowing explicit expression for the mean density of roots: 




47T A1A2 

with the coefficients 

K\ = —B cos 7 sin 7 + C cos 2 7 + D sin 2 7 
K2 = +B cos 7 sin 7 + C sin 2 7 + 13 cos 2 7 
K 3 = +B cos 2 7 + (C - D) sin 2 7 . 



(19) 



(20) 
(21) 
(22) 



Obviously, (p(r, ip)) can only be singular if at least one of 
the two eigenvalues Ai or A2 is zero. This condition leads 
to ab — c 2 . After some straightforward manipulations it 
can be written in the form Q(r, ip) = with 

N 

ip n - nip)r 



Q{r, ip) = C N C N sin 2 (^ m + m<p - 

n<m— 



2(n+m) 



(23) 



Thus, the points (r, ip) for which the average density of 
roots diverges are the zeros of the polynomial in Eq. (^3|) . 
However, Q{r, ip) is positive semi-definite. The only 
possibility of Q(rip) = is given by r = (which is 
always a solution and thus always a point of singular 
density), or by simultaneous vanishing of all coefficients: 
sm(ip m + mip — ip n — nip) — for all m, n. In the latter 
case Q(r, ip) will be zero for all r, implying immediately 
that lines of singular density can be only straight lines 
in the z-plane. On the other hand, assuming that ipk 
is a differentiable function of the index fc, one finds that 
<y5fc = —kip + const. with a fc-independent constant. Since 
the phases ip^ were chosen as constants, the only way to 
fulfill this equation is by ipk = ka + f3, ip = —a. For any 
other choice of the fc-dependence of the ipk, lines with 
more or less pronounced maxima of p(r, ip) may still ex- 
ist, but the singular character of the density is lost - 
with the exception of the origin. 

Above reasoning proves our claim that curves of sin- 
gular density are only possible if the phases ipk increase 
linearly with the index fc. This is exactly the case of 
the top [To, for which the symmetries T\ and T2 mani- 
fest themselves as singularities along straight lines on the 
complex plane, which correspond to great circles on the 
sphere. 

On the other hand, all the deviations from the above 
form result in a blurring of the sharp lines seen when 
plotting numerically obtained roots of random polyno- 
mials, irrespective of whether a particular symmetry of 
the possibly underlying physical system is still preserved 
or not. To demonstrate this effect we have analyzed 
random polynomials (|l|) with coefficients (Q) given by 
ipk = qk 2 /N . This assumption corresponds to the prob- 
lem induced by the generalized time-reversal symmetry 
(0) of the top 1/2- For q = (real coefficients a^) the 
distribution of zeros suffers a singularity along the real 
axis, while for larger value of q the clustering curve twists 
and acquires a finite width. This is visible in figures 2 
and 3 where we plotted on a complex plane zeros of 50 
random polynomials with N — 40 (part a) and the den- 
sity of zeros obtained according to Eq.([l9]) (part b). For 
q=0.2 the symmetry line already deviates from the real 
axis. For q = 0.5 a ridge in the density of zeros is still 
observed, at q ~ 1 the distribution of zeros is almost 
homogeneous. Interestingly, the qualitative character of 
the density does not change much with N. 

Let us mention here that the density of zeros of ran- 
dom polynomials (^) with fixed phases can be obtained 
using slightly different techniques proposed by Edelman 
and Kostlan [JhJ, Shepp and Vanderbei Jl5| or Prosen 
Moreover, the density of roots of some generalized 
random polynomials was recently discussed in p6[ . 
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in Essen, where this work has been initiated. This FIG - 3 - As in Fi § 2 for « = °' 5 - The concentration line of 

, , • ql ijuj-ucjiu the zeros is even more blurred than for q = 0.2. 
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FIG. 1. Mean entropy of eigenvectors compared with the 
entropy Hcue of the unitary ensemble drawn as a function 
of the perturbation parameter q for two models: "unitary" 
top Ui(A) and "orthogonal" top 1/2(0) with /i = 1.7, p = 6.0 
and spin length S — 40. The dashed line represents the value 
Hcoe/Hcue ~ 0.91. The inset shows the cumulative level 
spacing distribution P(s) obtained for both models out of 100 
operators U with fixed q = 2.0 and p varying from 6.0 to 12.0 
and compared to the Wigner surmises for both universality 
classes. 



FIG. 2. The distribution of roots of 50 random polynomi- 
als with quadratically increasing phases (q = 0.2) shown in 
part a) follows the analytically obtained density shown in the 
contour plot in part b). The concentration line of the zeros 
deviates from the real axis and is no longer a line of singular 
density. 
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